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ABSTRACT 

The N206 supernova remnant (SNR) in the Large Magellanic Cloud (LMC) 
has long been considered a prototypical "mixed morphology" SNR. Recent ob- 
servations, however, have added a new twist to this familiar plot: an elongated, 
radially-oriented radio feature seen in projection against the SNR face. Uti- 
lizing the high resolution and sensitivity available with the Hubble Space Tele- 
scope, Chandra, and XMM-Newton, we have obtained optical emission-line images 
and spatially resolved X-ray spectral maps for this intriguing SNR. Our findings 
present the SNR itself as a remnant in the mid to late stages of its evolution. 
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X-ray emission associated with the radio "linear feature" strongly suggests it to 
be a pulsar-wind nebula (PWN). A small X-ray knot is discovered at the outer 
tip of this feature. The feature's elongated morphology and the surrounding 
wedge-shaped X-ray enhancement strongly suggest a bow-shock PWN structure. 

Subject headings: supernova remnants — ISM: individual (SNR B0532-71.0) - 
X-rays: ISM — Magellanic Clouds 

1. Introduction 

The supernova remnant (SNR) B 0532—71.0 lies to the north and east of the H II region 
LHa 120-N206 1 in the Large Magellanic Cloud (LMC). The SNR was first identified in radio 
by Mathewson & Clarke (1973), who referred to it simply as N206. Although this is more 
properly the name of the larger H II complex, for simplicity we shall use the same notation 
in the following discussion. 

Ha images of N206 obtained from the Magellanic Cloud Emission-Line Survey (Smith et 
al. 1999) show largely shell-like emission, with complex filamentary structure along the limb. 
Williams et al. (1999) noted that these features were in sharp contrast to the X-ray emission 
as observed by ROSAT, which shows the X-ray emission brightening toward the center of the 
remnant. The striking difference between X-ray and optical morphologies led the authors 
to categorize N206 as a "centrally brightened" SNR, and to suggest that it might undergo 
the same physical processes as "mixed morphology" SNRs (Rho & Petre 1998), which have 
centrally brightened X-ray emission and a radio shell. 

Milne et al. (1980) studied N206 at 5 and 14.7 GHz, and found a spectral index a 
(as in S ~ u a ) over the remnant of about —0.33. This value is fairly flat for a SNR; it is 
on the borderline between the ranges for a pulsar wind nebula (PWN), for which spectral 
indices tend to be flatter than —0.4, and for a typical shell-type remnant, for which spectral 
indices tend to be between —0.8 and —0.3 (Trushkin et al. 2000). Low resolution (several- 
arcminute scales) and interference from the nearby H II region made features in the SNR 
interior difficult to examine. 

Of particular interest, therefore, are the observations by Klinger et al. (2002) of N206 
at 3 and 6 cm wavelengths (8.0 and 4.8 GHz), using the Australia Telescope Compact Array 
(ATCA) with resolutions of 1'.'8 and l'/l (HPBW), respectively. They saw emission over the 



1 Catalogue number from Henize (1956) 
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face of the entire SNR, somewhat brightened toward the limb. The spectral index calculated 
from their findings and those of other radio observations was — 0.20±0.07, in the normal 
range for filled-center SNRs. Klinger et al. (2002) also discerned a "peculiar linear feature" 
within the SNR: a narrow wedge of radio emission aligned in projection with the remnant's 
center. A spectral index map of the SNR showed this feature to have a similar index to the 
rest of the remnant, ~ —0.2. In the absence of high-resolution X-ray data, the authors could 
not isolate a point source, and therefore could only speculate that the feature was produced 
by "a low-mass star or compact object." Using the "opening angle" of the linear feature 
to calculate the Mach number of the presumed object's travel, and the length of the linear 
feature, the authors obtained an estimate for the SNR age of 23,000 yr. 



2. Observations 

2.1. X-ray Images and Spectroscopy 

X-ray observations of N206 were made with both the Chandra and XMM-Newton ob- 
servatories. Chandra observations used the Advanced CCD Imaging Spectrometer (ACIS), 
primarily the S3 back-illuminated chip. Datasets were sequence number 500327, observations 
3848 (33.1 ks) and 4421 (32.5 ks). Data were reduced following procedures recommended by 
the Chandra X-ray Center (CXC) and analyzed using the CXC's CIAO program and xspec. 
The two datasets were filtered for high-background times and poor event grades, resulting in 
a total "good time" interval of 65.6 ks. For each dataset, the 0'.'5 pixel randomizations were 
removed and the "subpixel resolution" technique 2 was applied, allowing us to more closely 
examine the structure on small spatial scales. 

The filtered datasets were then merged, and spectral results extracted from the merged 
file. As spectral analysis required us to choose fairly wide regions for reasonable counting 
statistics, the "subpixel resolution" files were not used for this purpose. A background region 
was taken from an annulus surrounding the SNR, and the spectrum of this background region 
was scaled and subtracted from the source spectra. Individual spectra for regions of interest 
selected from radio and X-ray images, and the corresponding primary and auxiliary response 
files, were extracted with the CIAO psextract script and analyzed in xspec. Spectra were 
rebinned by spectral energy to achieve a signal-to-noise ratio of 4 in each bin. 



2 This technique, developed by Koji Mori, allows users to improve the spatial resolution by 
~10% without decreasing the overall count rates. Documentation at http://cxc.harvard.edu/cont- 
soft /software / subpixel_resolut ion . 1 . 4 . ht ml 
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We received the pipeline-processed XMM-Newton data from the Science Operations 
Centre (SOC). Observations were made simultaneously with multiple XMM-Newton instru- 
ments; in this paper we will concentrate on the European Photon Imaging Camera (EPIC) 
MOS and pn detectors. The EPIC-MOS data were taken over two intervals, in 2001 Novem- 
ber (Observation ID 089210901, 41.4 ks) and 2001 December (0089210101, 14.5 ks). The 
latter observation also included an EPIC-pn exposure (12.0 ks). Initial reduction and anal- 
ysis were carried out using the Science Analysis Software (SAS) package provided by the 
SOC, with subsequent spectral analysis in xspec. 

The data were filtered to remove high background times or poor event grades, reducing 
the total "good" dataset time to 24.2 ks for the first observation (EPIC-MOS only) and 
12.0 ks for the second observation (EPIC-MOS and EPIC-pn). Images and spectra were 
then extracted from the filtered event files. Background regions immediately surrounding 
the SNR but free of point sources were used to produce background spectra, which were 
then scaled and subtracted from the source spectra. Joint spectral fits were performed with 
the data from the two EPIC-MOS observations. 



2.2. Optical Emission-line Images and Spectroscopy 

Images of N206 were obtained with the Hubble Space Telescope using the Wide Field 
Planetary Camera 2 (WFPC2) with the F656N (Ha), F673N ([S II] AA 6347, 6371 A) and 
F502N ([O III] AA 5007 A) filters. The southwestern side of the SNR was observed in Ha 
(3x800 s), [S II] (3x800 s) and [O III] (6x800 s), all at the same pointing and position angle, 
but the northeastern side of the SNR was observed only in Ha (3x800 s). The data were 
reduced using the iraf 3 program and STSDAS package. Multiple exposures were combined 
to remove cosmic-ray events, and the resulting files were bias-subtracted. The images were 
divided by exposure times to produce count-rate maps, which were then multiplied by the 
photflam parameter (provided in the image header) to convert these to flux-density maps. 
We used the synphot task to determine widths for each filter, and multiplied the flux- 
density maps by the filter widths to produce flux maps. The files for the WFC and PC were 
then mosaicked together for the final images. Images for [S II] and Ha were clipped at 3<r 
of sky background to reduce noise, and then divided to produce a calibrated [S II] /Ha ratio 
map. 

High-dispersion spectra of N206 were obtained with the echelle spectrograph on the 4 
m telescope at Cerro Tololo Inter-American Observatory (CTIO) from two observing runs, 



3 Image Reduction Analysis Facility, maintained by NOAO 
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2000 December 6 and 2004 January 14. Both runs used the 79 line mm -1 echelle grating 
in the single-order, long-slit observing configuration, in which a flat mirror replaced a cross 
disperser and a post-slit Ha filter (A c = 6563 A, AA = 75 A) was inserted to isolate a 
single order. The red long focus camera and the 2000x2000 SITe2K#6 CCD were used 
to record data. The 24 /zm pixel size corresponds to roughly 0.082 A along the dispersion 
and 0'.'26 perpendicular to the dispersion. Limited by the optics of the camera, the useful 
spatial coverage is ~3'. The spectral coverage is wide enough to include both the Ha line 
and the flanking [N II] AA6548, 6583 lines. A slit width of 250 fim (l'/64) was used and the 
resultant FWHM of the instrumental profile was 13.5 ± 0.5 km s _1 . The spectral dispersion 
was calibrated by a Th-Ar lamp exposure taken in the beginning of the night, but the 
absolute wavelength was calibrated against the geocoronal Ha line present in the nebular 
observations. The journal of echelle observations is given in Table 1. 



2.3. Radio Pulsar Observations 

We undertook new observations using the Parkes 64-m radio telescope at the Australia 
Telescope National Facility. We observed N206 for 12 hours on 2003 September 29 using 
the central beam of the 20-cm multibeam receiver at 1374 MHz on the Parkes telescope 
with the aim of detecting pulsed emission from a small X-ray source. Following the recent 
successful surveys for pulsars in radio nebulae at Parkes (Camilo et al. 2002a,b), we recorded 
total-power signals from 96 frequency channels which provide a bandwidth of 288 MHz for 
two polarizations every 1 ms. The sensitivity of the system was ~ 20/i Jy which is a factor of 
~2.5 more sensitive than an earlier, large-scale survey of the Large Magellanic Cloud (Fan 
2002). Off-line processing was carried out to search for any dispersed, periodic signal using 
standard procedures with the seek software. 4 No signal was detected. 



3. Physical Structure of the SNR 

The morphology of the SNR in the XMM-Newton and Chandra observations is com- 
plex, with particular structures revealing themselves variously under XMM-Newton's high 
sensitivity and Chandra's high resolution. In both cases (Fig. 1) diffuse emission can be seen 
over the entire face of the remnant, all the way out to the rim; and a significant increase in 
X-ray emission appears toward the center. In the XMM-Newton EPIC-pn and EPIC-MOS 
observations there is an extension of this central emission that corresponds well to the po- 



Source code and documentation are available at http://www.jb.man.ac.uk/~drl/seek. 
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sition of the radio "linear feature" . This apparent correlation is confirmed by the Chandra 
observations, in which X-ray emission is clearly seen along the length of the linear feature, 
broadening in the same region of the "wedge" as in radio (Fig. 2). 

An X-ray point source appears in the Chandra observation at 05 fe 32 m 03 s , — 71°00'51". 
No counterpart to this source is found in the optical or radio data. A search of the SIM- 
BAD database shows these coordinates to be within the error circle for Einstein source 2E 
0532.6—7102 (McDowell 1994). In the absence of features indicating an association between 
this source and the SNR, we presume it to be a background source; possibly, given its hard 
X-ray emission, an active galactic nucleus. We will not include this source in the subsequent 
discussion. 

The Hubble Space Telescope WFPC2 Ha mosaic of N206 in Figure 3 shows the circular, 
limb-brightened structure of this SNR. An intricate array of loop-like filaments extends from 
the outermost limb well within the remnant. Toward the southwest of the SNR the Ha 
emission becomes more prominent; this is unsurprising, as in this region the N206 SNR may 
overlap slightly with the larger N206 H II region. 

The [O III] emission, as seen in Figure 4c, tends to the outer edges of the filaments. 
[O III] is a useful "tracer" of shock fronts, due to the relatively narrow range of temperatures 
and densities in which this emission dominates. Thus the tendency of the [O III] to lie in the 
limbward direction of the filaments implies a fairly regular expansion. There do not appear 
to be any counterparts in Ha or [O III] to the linear feature seen in radio by Klinger et 
al. (2002). The Ha and [S II] morphologies, shown in Figure 4a and b, are quite similar, 
including the filamentary structure, as is typical for cooled (~10 4 K) post-shock gas. 

3.1. Hot Gas in the SNR's Interior 

Spectral fits to the X-ray data from the various instruments (Chandra ACIS, XMM- 
Newton EPIC-MOS and EPIC-pn) show that the emission from the SNR is dominated by 
thermal emission; a number of thermal emission lines are visible (Fig. 5). We therefore uti- 
lized thermal plasma models to determine the plasma parameters. As it is uncertain whether 
the plasma has reached ionization equilibrium, we fit both a non-equilibrium ionization model 
(NEI; "pshock" in XSPEC) and a collisional-ionization equilibrium model (CIE; "MeKaL" in 
XSPEC) 5 . We selected the following regions for further analysis: (Region 1) one comprising 



5 Details and references for the pshock and MeKaL models can be found at 
http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/node39.html and node40.html. 
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the entire SNR for comparison with other X-ray studies; (2) a region enclosing much of 
the SNR limb, which is expected to be dominated by emission from recently shocked gas; 
(3) a region enclosing the central X-ray brightening but excluding emission thought to be 
associated with the "linear feature" ; (4) a region covering an apparent wedge-shaped X-ray 
brightening north and south of the "linear feature", excluding a possible point source; (5) 
a region corresponding to the "linear feature" seen in radio; (6-8) three equal-area regions 
along this feature to examine possible changes over its length; (9) a small region surrounding 
a possible point source at one end of the "linear feature"; (10) a region covering emission 
immediately behind the possible point source, excluding a narrow area corresponding to the 
brightest radio emission from the "linear feature"; and (11) a region covering the aforemen- 
tioned radio-bright linear feature. A list of the regions used for X-ray analysis is given in 
Table 2 and shown in Figure 6. 

We chose two of these regions thought to include only thermal X-ray emission: "Outer 
Limb" (Region 2, Fig. 6b) and "Central" (Region 3, Fig. 6c). We compared the fits of simple 
thermal models to the data from each of the X-ray instruments used (Table 3). Best fits are 
obtained with sub-solar metal abundances, consistent with the ambient metal abundances 
in the LMC of about 30% solar (Russel & Dopita 1992). Fits to the absorption column 
density are reasonably consistent for regions with good statistics, giving N H = 3±1 x 10 21 
cm -2 . Fits to the temperature differ by a factor of two depending on whether CIE or NEI is 
assumed; CIE fits have kT = 0.23±0.01 keV and NEI fits have kT = 0.4±0.1 keV (Table 3). 
The fits to the ionization parameter r (= nt, where n is the hot gas density and t is the time 
since ionization) indicate that the hot gas deviates considerably from ionization equilibrium, 
with t = 3±1 x 10 11 cm~ 3 s. In general, the NEI model also provides a better statistical fit 
to the data. We will therefore refer primarily to the results of the NEI fits in the subsequent 
discussion. 



3.1.1. Overall SNR Characteristics 

In order to facilitate comparison of our X-ray results with other X-ray studies at lower 
resolution, we include the "Whole SNR" region (Region 1 in Table 2, Fig. 6a), which includes 
all X-ray emission from the remnant, including that from the area of the "linear feature". 
These fits (Table 4) largely fall within the range described above for the "Outer Limb" 
and "Central" regions, indicating that the SNR's emission is dominated by other sources of 
emission than that of the "linear feature." While an additional power-law component with a 
normalization of one-third that of the thermal plasma component slightly improves the fit, 
this improvement is not significant. The NEI fit yields an absorbed flux of 7±2 x 10~ 13 erg 
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cm 2 s 1 , an unabsorbed flux of 4±2 x 10 12 erg cm 2 s 1 , and a luminosity of 8±4 x 10 35 
erg s" 1 , all over the 0.3—8.0 keV range. 

The normalization of the NEI model fits to X-ray data can be used to calculate the rms 
electron density (n c ) of the hot gas within the SNR. We assume the remnant to be roughly 
spherical, with a radius of 21 pc, at a distance of 50 kpc. We further assume that the gas is 
composed of fully ionized hydrogen and helium (n e ~ 1.2nn). Then, using the normalization 
from NEI fits to the Chandra ACIS data for the "Whole SNR" region (Table 4) we find a 
density in the hot gas of n e = 0.24 ± 0.05 / llo ]/ 2 cm~ 3 . Using the "Outer" and "Central" 
regions (Table 5-6), which do not include emission from the area around the "linear feature," 
for this calculation gives a very similar density of n e = 0.23±0.09 / h J/ 2 cm -3 . The parameter 
/hot is a volume filling factor for the hot gas. If we assume /hot=l, to reflect the centrally 
filled X-ray morphology of the remnant, the hot gas density is n c = 0.24 ± 0.09 cm -3 ; if the 
gas only partially fills this volume, the density will rise in inverse proportion to the square 
root of this filling factor. 

From this derived density of the X-ray emitting material, we obtain a mass of 5.3±0.2 x 
10 35 g> or 270±10 fl^t M©. Presuming that the hot gas fills the entire remnant (/hot = l)> 
this mass is greater than that expected from SN ejecta alone, emphasizing that N206 is an 
older SNR whose emission is dominated by swept-up material from the surroundings; and 
supports the LMC-like abundances found from the spectral fit. Using the temperature and 
density from spectral fits to the X-ray data, we find a thermal pressure, P t h = nkT, of 
3.3±0.4 x 10~ 10 f h ot 2 dyne cm -2 in the hot gas. As the analysis in §3.2 will show, this 
is significantly greater than the thermal pressure we find for the cool shell. We can also 
calculate the thermal energy in this gas, E t b. = 1.5nfVkT. This gives a thermal energy 
of 6±1 x 10 50 erg for the hot gas. Finally, we can estimate the shock velocity under 
the (somewhat dubious) assumption that the bulk of the X-ray emission is being generated 
at the current shock front, using kT = (3/16)//i> 2 shock1 with the reduced mass /i = 0.61 (for 
He:H = 1:10). If we presume the newly shocked gas to be at 0.45 keV, we find that such a 
temperature would result from a shock speed of 620 km s -1 , implying an overall expansion 
velocity of 470 km s -1 . For simplicity, we have assumed in the above calculations that the 
ions have equilibrated with the electrons; although this assumption may well not be valid, 
it provides us with at least initial estimates. The derived energy and pressure are similar to 
those for other SNRs in the adiabatic stage of evolution. 
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3.1.2. Specific Hot Gas Features 

N206 is clearly more complex than the simple shell approximated above. It features an 
increase in surface brightness toward the remnant center, a "wedge" of emission appearing 
to surround the "linear feature", and the feature itself, terminating in a possible compact 
source. The properties of these features will be discussed in more detail in §4. 

The central brightening on the SNR is quite puzzling. The spectrum of this region 
(Table 6) appears similar to that for the limb of the remnant, with slightly higher abundances, 
particularly oxygen. There is no evidence for a strong power-law component in this region; 
fits with such a component require it to be normalized to a small fraction of the thermal 
component's emissivity. 

Toward the eastern side of the remnant, surrounding the "linear feature," the X-rays 
show a slightly brighter broad "wedge" of emission (Region 4 in Table 2, Fig. 6c). Fits 
to the spectrum of this region (Table 7; this spectrum does not include contributions near 
the "compact source" ) suggest a markedly higher plasma temperature than elsewhere in the 
SNR. While a nonthermal contribution might artificially skew a thermal model fit toward 
higher temperatures, fits which included nonthermal emission over a tenth of the thermal 
emissivity could be statistically ruled out at the 90% level. 

We also examined a region covering the area of the "linear feature" , as seen from radio 
images (Region 5 in Table 2, Fig. 6d). We defined a hardness ratio (H— S/H+S) such that 
S=0. 3-1.0 keV and H=l. 0-8.0 keV (Table 2). Using these hardness ratios to compare this 
region to others throughout the SNR, we find that, while the contribution of soft-X-rays still 
dominates, the proportion of hard X-rays in this region is greater than elsewhere in the SNR. 
This increase in the hardness ratio is consistent with the combination of thermal emission 
from the SNR and emission from a harder source, as for example the small-diameter bright 
source at the tip of this region. While a localized increase in temperature could also explain 
this hardening, it would be quite coincidental that this increase appears only in this portion 
of the SNR that also shows the small-diameter X-ray source and elongated radio feature 
(Fig. lc-e). 

We performed spectral fits to data from Region 5, though we note that these fits are 
somewhat limited by the relatively low number of counts. We find that either a thermal 
plasma model or a power-law model can provide a statistically adequate fit to this source 
(Table 8). The "best" fit arises from a combination of nonequilibrium thermal plasma and 
power-law models, providing roughly equal contributions to the spectrum. Although the 
improvement in the fit is not statistically significant, such a combination is consistent with 
the scenario inferred from the hardness ratios. To examine the detailed properties of the 
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"linear feature" , we used the Chandra ACIS data to analyze emission from smaller segments 
along that feature (Regions 6-8 in Table 2, Fig. 6e). These fits, summarized in Table 9, 
should be considered as preliminary estimates only. 

3.2. The Cool Shell of the SNR 

Figure 4d shows an [S II] /Ha ratio map of the southwestern side of N206. [S II] /Ha 
ratios range from 0.7 to 1.2 across the SNR, typical for shocked gas, with an average value 
of 0.9. This is initially somewhat surprising, as it indicates that there has not been much 
ionization of material from the nearby OB association. However, H I maps of the vicinity 
show an H I shell, with the N206 H II region largely situated within the central cavity, so 
that the H I shell walls probably block most of the ionizing radiation before it can reach the 
SNR (Dunne, Chu & Stavely-Smith 2005). 

We measured the average Ha surface brightness of filaments in the flux-calibrated 
WFPC2 image of N206 to be (0.8-1.8) x 10" 17 erg cm" 2 s" 1 pix -1 , or (0.8-1.8) x 10" 15 erg 
s" 1 cm" 2 arcsec" 2 . Assuming a cool shell temperature of 10 4 K, this surface brightness im- 
plies an optical emission measure of 680± 260 pc cm -6 . We presumed the average filament 
thickness along the line of sight (20" — 30") to be equal to its width perpendicular to the 
line of sight, and used this as a representative number for the path length £ through the 
warm ionized gas. Using this C, we calculate an rms electron density in the shell of about 
10±4 cm" 3 . If we presume the SNR to be in the adiabatic "point-blast" stage of expansion 
(Sedov 1959), and the cool shell to be representative of the ambient ISM, we would expect 
the ambient ISM density to be roughly one-quarter of that in the shell, or about 3±1 cm -3 . 

The optical echelle data (Fig. 7) show several lines in the Ha spectral region, includ- 
ing the narrow geocoronal Ha (6562.85 A) and telluric OH 6-1 P2(3.5) 6568.779 and 6-1 
Pl(4.5)e/f 6577.183/. 386 lines (Osterbrock et al. 1996), as well as broader lines correspond- 
ing to Doppler-shifted nebular Ha emission toward the SNR. The latter include both a 
velocity component constant along the slit, showing the systemic velocity of the SNR, and 
the characteristic bow-shaped pattern deviating from this systemic velocity, showing motions 
within the expanding gas. 

In order to measure the systemic velocity, we extracted a velocity profile from a region 
13" wide outside of the emission of the SNR expansion pattern. The profile showed two 
components, one of which was identified as the telluric OH line. The other component is 
the Ha line with a measured Doppler shift of 241 ±4 km s -1 ; as the SNR expansion pattern 
appears to converge to this component, we take this as the systemic velocity for the N206 
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SNR. For comparison, we note that the nearest position to N206 in the H I maps of Rohlfs 
et al. (1984) shows velocity components at 240±7 and 262±26 km s -1 . 

To characterize the expansion of the N206 SNR, we measure the velocity offsets (Av) 
from the systemic velocity in both the blue and red directions. As the material at the 
forefront of the SNR expansion may be reasonably expected to show the highest velocity, 
we take the greatest of these velocity offsets to represent the overall expansion velocity of 
the remnant. The greatest offset in the blue direction (Atonic) * s — 202±5 km s" 1 , while the 
largest At> re d= +193±5 km s _1 . We therefore estimate the expansion velocity f CX p = 202±5 
km s _1 . Using this value for expansion, and assuming that the SNR is in the Sedov phase, 
we would expect the shock velocity f s hock = 4/3 t> cxp = 270±5 km s" 1 . It should be noted 
that the errors given here are only the random errors of the measurements. 

It is possible that some of the optically emitting material is too faint for detection 
in these observations. If the highest-velocity material is undetected, the actual expansion 
velocity may be higher than our estimate. For example, Chu & Kennicutt (1988) found 
a higher expansion rate of 250 km s _1 ; the discrepancy may be due to the difficulty of 
discerning motions in the faint outer material. To an extent, this may also apply to the large 
discrepancy between the shock velocity calculated from the temperature of the hot gas in 
§3.1.1, and that measured from the warm ionized gas. More likely, however, this discrepancy 
is due to two factors. (1) The X-ray temperature may not be a reliable representation of 
the newly shocked gas, particularly given its large deviation from a limb-brightened shell 
morphology. (2) We expect the X-ray emission to arise from areas where the shock front 
is moving through diffuse gas, while the optical emission is expected to be generated where 
the shock is moving through higher density gas, as seen in the relative densities of the hot 
gas to the warm ionized gas. Thus, we expect that the shock within the optically-emitting 
clumps has been somewhat slowed by its progress through this denser material. The overall 
expansion velocity from the blast wave is probably intermediate between the measured optical 
expansion of 202 km s _1 and the calculated expansion from X-ray temperatures of 470 km 
s" 1 . 

Using the shell electron density of n=10 cm" 3 calculated above, and presuming nu e /nu = 
0.1 and singly ionized He (n e = I.Iuh), we can infer the mass of gas in the shell according 
to M — nmnUsheii, where Usheii is the shell volume. From the filamentary structure at the 
edge of the shell, we measure an average shell thickness of about l''l±0 / /5 (0.27±0.13 pc at 
a distance of ~50 kpc to the LMC). If we assume a simple spherical shell of this thickness, 
we calculate a shell mass of 9±6 x 10 35 g, or 460±300 M Q . This is almost certainly an 
underestimate, as the filamentary structure shows a much greater extent of cool gas than 
such a simplistic scenario. If instead we presume the cool gas occupies a volume filling factor 
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of as much as 0.1, we obtain a mass estimate of 2.4±1.7 x 10 36 g, or 1200±800 M Q . 

Using these as low- and high-end estimates of the mass range, and the expansion velocity 
above, we calculate the kinetic energy in the cool shell, 0.5M shcll t^ xp , to find values of 2±1 x 
10 50 erg and 5±3 x 10 50 erg, respectively. One can also use the density calculated above to 
calculate the thermal pressure in the shell, P = nkT, where k is the Boltzmann constant 
and T is the temperature, presumed here to be 10 4 K. This equation gives a pressure of 
3±1 x 10~ n dyne cm -2 . Again, these values for kinetic energy and shell pressure are typical 
for middle-aged SNRs. The physical characteristics of the N206 SNR are summarized in 
Table 10. 



4. Specific Radio and X-ray Features 

The high spatial resolution of Chandra's ACIS allows us to examine specific regions in 
more detail. We have isolated three features that appear to be of particular interest: X-ray 
emission associated with the radio- identified "linear feature;" the X-ray knot and surrounding 
emission which coincides with the very tip of that linear feature; and the region of enhanced 
X-ray surface brightness toward the center of the SNR. Below, we describe each of these 
features individually and then discuss the possible origins for these features. 

4.1. Central Emission 

The interior X-ray emission is brightest from an irregular region extending N— S for half 
the diameter of the remnant. This central emission is clearly dominated by thermal X-rays, 
with prominent line features including the blends of He-like lines of Mg and Si. The emission 
is largely soft (< 2 keV) and thermal plasma fits to data from this region give parameters 
for Ah and kT similar to those found for the SNR as a whole. However, there does appear 
to be some difference in abundance distributions between the central area and the outer 
regions of the SNR (excluding emission associated with the "linear feature"). The oxygen 
lines in the central region are more prominent, with respect to the iron blends (Fe L), than 
in other regions of the SNR (Table 6). These differences should be treated with caution, 
as there are large uncertainties in the abundance determinations; however, the difference in 
oxygen abundance between the two regions within the SNR is significant when compared to 
the 90% error ranges for this quantity. If the central emission is due to "fossil radiation," 
i.e. gas that was shock-heated during earlier phases of the SNR expansion, we would expect 
that this gas would have a higher proportion of ejecta to swept-up matter than the more 
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recently-shocked gas, and therefore that the higher oxygen abundance reflects the presence 
of oxygen-rich ejecta, as seen in the much younger mixed-morphology SNR 0103-72 in the 
SMC (Park et al. 2003). The presence of O-rich material cannot be confirmed from Hubble 
Space Telescope [O III] maps. This probably reflects the fact that the X-ray emission which 
shows the oxygen excess is primarily from the hot central cavity, where the temperatures are 
too high and the densities too low for [O III] to show up. 

It is worth noting that in addition to the similar oxygen-rich and mixed-morphology 
SNR 0103-72 (Park et al. 2003), the oxygen-rich SNR N132D shows, likewise, a higher ratio 
of oxygen (0 7+ ) to other elements in the remnant interior than at the bright rim (Behar et 
al. 2001). Nucleosynthesis models indicate that a high oxygen abundance within a SNR's 
ejecta is the result of a Type II SN (e.g., Tsujimoto et al 1995). The examples of "mixed- 
morphology" SNRs which show signs of enhanced oxygen abundances, therefore, may imply 
that the "mixed-morphology" SNRs are likely to have originated from Type II SNe. 

As with other mixed- morphology SNRs (Rho & Petre 1998), the expanding shock of 
N206 has slowed to the point where bright X-rays from the limb do not dominate the overall 
emission from the remnant. Thus, it is reasonable to think that N206 follows the typical 
pattern for mixed-morphology SNRs, wherein the central emission is dominated by fossil 
radiation from earlier large-scale shock heating during the SNR's evolution. 

4.2. Linear Feature 

An elliptical region (65" E— W x 21" N— S) surrounding the "linear feature" in radio 
was selected from the 6 cm radio image (Klinger et al. 2002), and was used to extract spectra 
from the corresponding X-ray data (Region 5). Comparison of the hardness ratio for this 
region to those elsewhere in the SNR, as discussed in §3.1.2, suggests the presence of a 
hard component in addition to the soft emission seen throughout N206. A combined NEI 
and power-law model with kT = 0.4 ± 0.2 and T = 2.2 ± 0.2 provided the best fit. Fits 
to these spectra are summarized in Table 8. (Note that the value of r obtained for this 
fit would indicate a plasma in collisional ionization equilibrium; we retain the NEI model, 
however, for consistency in intercomparison with the other fits.) Using these fitted values, 
we obtained an absorbed X-ray flux of 4±lx 10~ 14 erg cm -2 s _1 over the 0.3-8.0 keV range, 
which corresponds to an unabsorbed flux of lx 10~ 13 erg cur 2 s -1 . At the LMC distance 
of 50 kpc, this gives a luminosity of 3x 10 34 erg s -1 over this energy range. 

To further examine the nature of this linear feature, a series of three equal-area regions 
(20" x 15") E— W along the feature were identified in Chandra images for further spectral 
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examination (Regions 6-8 in Table 2). Comparison of the hardness ratios defined above for 
regions 6-8 shows a pronounced shift from regions dominated by high-energy photons to 
those dominated by low-energy photons. 

We also performed spectral fits for these regions (Table 9). It should be stressed that, 
due to the low numbers of counts, these spectral fits are quite uncertain. Cited values of 
Xred are l° w due to the large error bars for the spectral bins; however, further binning was 
deemed undesirable due to the loss of remaining spectral information. However, as with the 
hardness ratios discussed above, the fits illustrate the shift from harder to softer emission 
along the feature. 

In order to see whether the relative fluxes of thermal to nonthermal emission would 
change substantially as one moved from the remnant center toward the X-ray bright knot 
at the east end of the linear feature, we adopted a dual approach to these spectral fits. 
Initially, we fit a two-component (power-law and NEI) model jointly to all three regions, and 
determined the flux in each component for each region. As an alternate approach, we fit the 
power-law component to the high-energy end of the spectrum only (>2 keV) and fixed those 
parameters, then fit the joint model to the low-energy end of the spectrum. Again, the flux 
in each component was determined. The results, summarized in Table 11, indicate that the 
X-ray emission along the linear feature is more strongly dominated by nonthermal emission 
the further one moves eastward (closer to the knot) along the feature, and dies off on the 
west side. While this result is largely qualitative, given the errors, we argue that the increase 
in "hard" emission toward the knot is most likely to result from a nonthermal source. 

The nonthermal spectra of the X-ray emission and radio emission in the "linear feature" 
imply the presence of a pulsar-wind nebula (PWN). The most likely source for nonthermal 
X-rays within a SNR is synchrotron radiation. While nonthermal X-ray emission may occa- 
sionally be generated by a fast SNR shock at the outer limb of a SNR, both the relatively 
slow expansion of this SNR and the shape of the emitting region of the nonthermal emission, 
i.e. over a small spatial segment perpendicular to the shell, argue against this scenario. 
Likewise, it seems unlikely that the motion of a compact source would be sufficient to pro- 
duce a bow shock capable of generating significant synchrotron X-rays. The most plausible 
explanation, then, is that the acceleration of particles by a "hidden" pulsar is the source of 
the nonthermal emission (in both radio and X-rays) seen in the linear feature. 

The westward decrease in the ratio of nonthermal to thermal X-ray flux along the linear 
feature suggests that it is a PWN generated by an eastward moving pulsar. The lifetime 
of the high-energy particles generated in a PWN is comparatively short; the radio-emitting 
particles persist and so trail farther behind the X-ray extent. The radio pulsar observations 
were unable to find any periodicity greater than 2 ms to the radio signal from the area of 
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the "linear feature" . It is by no means uncommon for a PWN to be detected in radio and 
X-rays without a radio point-source counterpart; for example, N157B, also in the LMC, was 
deduced to have a PWN long before the X-ray pulsar was discovered, and a radio counterpart 
for that pulsar has yet to be found (Wang et al. 2001). 

Another possibility that must be considered, of course, is that these features are as- 
sociated with a background source rather than the SNR. Klinger et al. (2002) discuss this 
scenario, but conclude that this is unlikely. The presence of extended nonthermal X-ray 
emission provides an additional reason to believe the feature is not a background source, as 
such a feature would be highly unusual in a background galaxy. 

If we accept the premise that the "linear feature" is a PWN, we may anticipate that 
in part, the morphology of the PWN is created by a bow-shock structure. The dynamics 
of a rapidly moving pulsar and evolving SNR lead to the prediction of the formation of a 
bow shock when the SNR expansion speed has decelerated sufficiently for the pulsar's motion 
relative to the shocked SNR material to become supersonic (e.g. van der Swaluw et al. 1998). 
At this point the PWN is deformed, leading to a substantial offset between the pulsar and 
the center of the PWN; this is certainly consistent with the observed bright knot at the outer 
tip of the "linear feature" presumptive PWN. 

Numerical simulations predict that for a remnant in the Sedov stage, this should occur 
when the distance traveled by the pulsar from the center of the SNR, -Rpsr, approaches the 
SNR radius -Rsnr according to -Rpsr/-Rsnr ^ 0.677 (van der Swaluw 2004; van der Swaluw 
et al. 1998). Presuming that the geometric center of the SNR shell accurately represents 
the site of the SN, the current position of the small-diameter X-ray source is such that the 
transverse component of -Rpsr is 0.85 -Rsnr, well over the point for bow shock formation. 
If the direction of pulsar motion is not perpendicular to the line of sight, -Rpsr increases 
accordingly. These findings are in general accord with those of Klinger et al. (2002), and 
their assumptions in their estimation of the pulsar velocity and age. 

The fact that emission from the "linear feature" in X-ray and radio does not dominate 
the overall emission from the remnant presents an additional question. Its spectral index 
would place it among the most flat radio SNR shells, although overlapping (within the error 
bars) with remnants such as IC443 (Kawasaki et al. 2002). It is difficult to see why the 
SNR as a whole should show a radio index typical of filled-center remnants (Klinger et al. 
2002) if the putative PWN is largely confined to this feature. However, the radio emission 
from N206 is faint compared to that from other SNRs, and therefore the uncertainties in the 
spectral index determination are considerable. In addition, residual emission from the N206 
H II region may contaminate the radio emission, adding additional uncertainty. 
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4.3. Evidence for a Compact Source 

In the Chandra images, a distinct knot of X-ray emission appears to the far eastward 
end of the linear feature, near the SNR limb. To investigate whether this knot is consistent 
with an unresolved source, we generated a point source of comparable X-ray flux using the 
MARX simulator, and compared the profile of this simulated source to that of the knot. We 
conclude that the knot is a small diffuse region of dimension 2". This can also by seen by 
comparison with the appearance of the moderately bright point source ~50" SW of the knot 
(identified above with 2E 0532.6-7102). If a point source is embedded within the knot, only 
<50% of the counts from the knot could be due to that point source. 

We extracted counts within a 2" radius of the knot from the merged Chandra ACIS 
events file (Region 9 in Table 2, Fig. 6e), obtaining a count rate after background subtraction 
of 7.9 x 10~ 4 ct s^ 1 , or 52 counts over the 65.6 ks exposure time, an insufficient quantity to 
obtain a meaningful spectrum. Fig. 1, however, shows that the knot shows noticeably harder 
emission than that from the surrounding SNR. Once again using the hardness ratio defined 
above, (Table 2), we find a hardness ratio of 0.5 for the knot, in contrast to hardness ratios 
of -0.45 and -0.58 for the "outer" (Region 2; SNR limb excluding the putative PWN) and 
"central" (Region 3) areas. 

To obtain a first-order estimate of spectral properties for this knot, we fixed Ah (and 
abundances, for the thermal case) to the values determined by fits to the X-ray emission from 
the rest of the SNR. The data are consistent with a thermal plasma model with temperatures 
> 1 keV, or a power-law model with a spectral index Y between 2 and 3. Using a power- law 
model with V = 2.4, we obtain an estimate for the absorbed flux of 6x 10 -15 erg cm -2 
s _1 , an unabsorbed flux of lx 10" 14 erg cm -2 s -1 , and a luminosity of 3x 10 33 erg s" 1 at 
the LMC distance, all over the 0.3-8.0 keV energy range. The maximum luminosity for an 
embedded point source is therefore ~1.5x 10 33 erg s" 1 over this range. 

An elliptical region (5"x8") immediately west of the knot, along the linear feature, 
yielded a count rate of 1.4 x 10~ 3 ct s _1 . This region has a very similar spectrum to that 
from the knot, suggesting that the spectrum determined from the knot is dominated by 
emission from its immediate surroundings. Combining the spectra from these two regions 
allowed us to increase the signal-to-noise ratio for a slightly better spectral fit. While these 
fits would still not allow us to rule out a thermal plasma interpretation at 90% confidence, 
they confirm that a relatively high temperature kT > 1 keV, is required. More plausible 
is a power-law fit, with T = 2 ± 1, roughly consistent with that measured from the "linear 
feature" as a whole. While fits to smaller regions along the linear feature suggest a possible 
spatial variation in the power-law spectral index, the sensitivity of the X-ray observations is 
insufficient to make this determination at a statistically significant level. 
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A brief analysis of the X-ray power spectrum shows no evidence for periodic emission 
from an embedded pulsar. However, the timing resolution of these observations, 3.2 seconds, 
would not be sufficient to detect the pulsations from a typical pulsar. In addition, the 
estimated maximum of 25 counts from a point source is insufficient to determine a period 
for its emission, even if the timing resolution were available. 



Proceeding westward from the knot, X-ray emission associated with the "linear feature" 
broadens into a wedge of brighter X-ray emission which then appears to merge with the 
central emission. Faint emission connects this patch to the compact knot close to the eastern 
rim. We interpret this connection as a "bow shock," clear in the south and barely discernable 
in the north. This bow shock merges with the radio "linear feature" close to the knot. A 
bow shock in front of a hypothetical moving pulsar might re-heat the material through which 
it moves, leaving a trail of "fossil radiation" back to the remnant center and merging there 
with fossil radiation from much earlier fast-moving shocks. 

One approach to these features would be to consider much of the X-ray emission between 
the bright linear feature and the central region as the actual boundary of the bow-shocked 
gas. In this case we would have a much larger opening angle for this structure than that 
found by Klinger et al. (2002). Using the X-ray opening angle of about 57° ± 5°(1.0 ±0.1 
rad) we obtain a Mach number Ai = sin(#/2) -1 = 2.1 ±0.5. This is considerably lower than 
than the M. = 9 value of Klinger et al. (2002), and closer to the theoretical value for a pulsar 
moving through a SNR interior (M. 3.1, van der Swaluw et al. 1998). 

We can calculate the sound speed of the hot gas within the SNR according to 



where 7 is the adiabatic index, here taken as 5/3; k is Boltzmann's constant; T e is the electron 
temperature; \i is the reduced mass; and m H is the hydrogen particle mass. Presuming the 
reduced mass ji = 0.61, and using the temperature found above, we obtain a sound speed 
in the hot gas of 230 ± 110 km s _1 . According to M = u/c s , where u is the motion of 
the pulsar relative to the SNR material, this gives a mean relative motion of w=480 ± 240 
km s -1 . 



4.4. Possible Bow Shock 




If the pulsar is, as it appears, moving radially outward from the SNR center, we can 
presume its motion to be parallel to that of the expanding SNR material around the pulsar. 
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In order to obtain a mean relative velocity of w=480 km s _1 with material moving at the 
expansion velocity t> cxp = 202 km s -1 (slow-expansion case), we would require the pulsar to 
be moving at a speed of roughly 680 km s -1 with respect to the SNR center. If instead we 
use the X-ray derived expansion speed of t> exp = 470 km s" 1 (fast-expansion case), the pulsar 
must be moving with a speed of ~950 km s _1 with respect to the SNR center. If the motion 
is entirely perpendicular to the line of sight, and the pulsar began at the geometric center 
of the SNR, it would have taken the pulsar about 28,000 yr to reach its current transverse 
distance of 19.4 pc in the slow-expansion case, and 20,000 yr in the fast-expansion case. 
These can, of course, be affected by the viewing angle; but if we assume the pulsar to be 
interacting with the SNR interior, its maximum distance of travel is the SNR radius of 21 
pc (at a 67° angle to the line of sight), which puts its time of travel at 30,000 yr in the 
slow-expansion case and 22,000 yr in the fast-expansion case. 

In our suggested picture, the bow-shocked gas is primarily thermal, representing material 
within the SNR that is shocked to higher temperatures by the encounter with the bow shock. 
Presumably, the nonthermal X-ray emission is associated with the elongated PWN, as would 
be the radio emission. Thus we would expect a broad outer "cone" of thermal (shocked SNR 
material) emission, with an interior "cone" of nonthermal (shocked PWN) emission. Such 
a scenario is in accord with hydrodynamic models of bow-shock PWN, in which "the Mach 
cone should manifest itself only in the outer bow shock" (Gaensler et al. 2004, and references 
therein). The shocked PWN material forms a trail of synchrotron emitting particles opposite 
the direction of pulsar motion, creating a "cometary" or "linear" morphology. 

To examine this picture for consistency with the data, we fit a simple thermal plasma 
plus power-law model combination to data from two regions within the linear feature: the 
thin "trail" region, corresponding to the narrow location of the brightest radio feature (Re- 
gion 11 in Table 2, Fig. 6f), and the "bow" region, comprised of emission around the linear 
feature excluding the aforementioned "trail" (Region 10, Fig. 6f). Both regions excluded 
emission within 6" of the presumed pulsar. The best-fit model combination gave thermal 
and power-law parameters similar to those given in previous sections (Nh ~ 2.6 x 10 21 cm -2 , 
T rsj 0.26 keV, T ~2.3). Calculating the fraction of the total flux for each model compo- 
nent, we find that the "trail" region shows 96% nonthermal flux and only 4% due to thermal 
plasma, while the "bow" region shows 49% thermal flux and 51% nonthermal flux. Errors on 
these flux ratios, based on the differences in flux due to uncertainties in the fitted parameters 
at the 90% confidence level, are of order 10%. 

Alternately, it is possible that some of the X-ray brightening we consider to be associated 
with the radio "linear feature" (excluding the X-ray knot) is in fact only a surface brightness 
fluctuation, and does not accurately represent the "Mach cone" of the presumptive moving 
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pulsar. If the "Mach cone" is more tightly confined than our X-ray based estimate, the 
resulting Mach number may be somewhat higher, up to a value of about Ai = 4 when 
one considers only the angular extent of the bright X-ray emission immediately surrounding 
the X-ray knot. While it is also possible that the emission from these features is due to 
the PWN alone, with no substantial bow-shock emission, the highly elongated morphology 
and off-center pulsar position argue strongly for the presence of a significant bow-shock 
contribution. 



5. SNR age 

The ages of supernova remnants are notoriously ill-characterized. We consider a number 
of different approaches to this age, in order to find the range of reasonable age estimates. 
From the fitted X-ray parameters, several estimates are possible. Combining the derived 
hot gas density with the fitted ionization parameter r gives an upper limit age estimate of 
~40,000 years for this remnant. However, an age this large would be difficult to reconcile with 
the fact that the hot gas in the SNR has apparently not yet reached ionization equilibrium. 
Lowering the filling factor shortens this age estimate significantly. A filling factor of 0.25, 
for instance, gives an age estimate of 23,000 yr, similar to the estimate based on the radio 
"linear feature" by Klinger et al. (2002). 

We can also make another estimate from the Sedov relations, in which the relative 
fractions of thermal and kinetic energies to total energy are constant. Taking the theoretical 
relation of E t ^ 0.7E , our estimate for thermal energy (for /h ot = 1) would yield an initial 
explosion energy E of 8 x 10 50 erg. Assuming the density behind the shock is a factor of 4 
greater than that of the unshocked ISM, we obtain an ISM mass density po of 1.2 xl0~ 25 g 
cm -3 . Then from the Sedov relation R = 1.17(i?o/po) 1 ^ 2/ ' 5 > we have an age of only 9,000 yr. 
If the X-ray filling factor is lower than 1, the age estimate will rise accordingly; / ho t = 0.25, 
for instance, gives an age of 17,000 yr. If we instead presume the overall ambient density 
to be represented by the current filaments of warm ionized gas, we obtain an ISM mass 
density of 5.3 xl0~ 24 g cm -3 . The Sedov relation then gives an age as high as 57,000 yr, 
demonstrating the critical role played by the ambient density assumption. 

Using the the simple expansion relation t = i]R/v cxp , and assuming Sedov expansion 
(?7=0.4), we can calculate the age of the remnant based on the expansion velocity. The 
expansion velocity of 470 km s _1 derived from the X-ray temperature gives a remnant age of 
17,000 yr. The expansion velocity obtained from optical echelle spectroscopy of 202 km s _1 , 
however, leads to a calculated age of 41,000 yr. Again, however, this must be considered an 
upper limit, as this value of t> cxp is, as stated above, a lower limit to the blast- wave expansion. 
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We can try to narrow down these estimates by considering the consistency of the entire 
picture. Ages above ~40,000 yr are clearly ruled out by the ionization timescale of the X-ray 
gas and by the observed optical expansion. The value of 9000 yr from energetics arguments 
is highly dependent on the assumption of ambient density; a better lower limit might be 
that of 17,000 yr from the expansion velocity derived from the X-ray temperature. When we 
consider the estimates of pulsar motion (§4.4) we find further limitations: the minimum time 
for the pulsar to reach its current location is ~20,000 yr, while the maximum is ~30,000 yr. 

If we choose a velocity intermediate between the X-ray and optically derived values (e.g., 
~300 km s _1 ) and a hot gas filling factor of 0.25, we can reconcile most of these numbers. 
The ionization timescale for the hot gas then gives an approximate age of 23,000 yr. Pulsar 
travel times give ages (depending on viewing angle) between 24,000 - 26,000 yr. Simple 
Sedov expansion at this velocity gives an age of 27,000 yr. While arguments from Sedov 
energetics would, for this same filling factor, give us a somewhat lower age of 17,000 yr, this 
estimate is highly dependent on the assumption used for ambient density, a very uncertain 
parameter. We therefore suggest the SNR's age most likely falls in the range between 23,000 
- 27,000 yr, with "hard" limits of 17,000 yr - 40,000 yr. 



6. Integrated Picture 

We have observed the SNR in N206 with high-resolution optical and X-ray instruments, 
and analyzed the results in concert with additional radio data and optical echelle spectra. 
We find it highly probable that N206 is the result of a Type II SN, due to its proximity to 
other massive-star phenomena, enhancement of oxygen abundances, and the presence of a 
probable compact object. 

We use these data to calculate overall properties for the physical components of the 
SNR, and find these typical of a middle-aged SNR in the adiabatic stage. The remnant is 
over-pressured, and the bulk of the energy budget still resides in thermal energy from the 
hot interior. Evidently, therefore, the hot gas within the SNR still plays a significant role in 
the remnant's development. 

This remnant is particularly unusual for the very different characteristics it displays in 
different wavelength regimes. The optical morphology is limb-brightened and highly filamen- 
tary; the radio morphology is center-filled, with diffuse emission over the remnant's face, but 
some SNR shell structure; and the X-ray morphology appears somewhat centrally bright- 
ened. The best categorization for this remnant appears to be that of "mixed-morphology" 
(Rho & Petre 1998), but the picture is complicated by the presence of radio and X-ray 



- 21 - 



emission near the "linear feature" . 

We analyze the X-ray data for the area surrounding the "linear feature" seen in radio. 
We find a small, hard X-ray source located at the tip of the radio feature, with a surface 
brightness profile consistent with the presence of an embedded compact source. Emission 
from this source and its surroundings is nonthermal, with a power-law index similar to that 
seen in Crab-type objects. The ratio of nonthermal to thermal X-ray flux decreases with 
increasing distance from this source. We conclude that the most probable scenario for this 
feature is a pulsar moving at moderate velocity through the surrounding SNR. This creates 
a bow-shock structure in the direction of motion, deforms the surrounding PWN, and leaves 
behind a trail of synchrotron emission along the line of travel. 

However, it would be difficult to attribute most, or even the majority, of the central 
X-ray emission to that linear feature. The X-ray emission from the central region is clearly 
dominated by thermal line emission, rather than nonthermal continuum. The lifetime of 
relativistic particles energetic enough to generate X-ray emission is not sufficient for those 
particles to be responsible for the central X-rays, as indicated by the fact that the radio emis- 
sion, generated by longer-lived particles, shows no particular increase in this central region. 
This indicates, somewhat surprisingly, that N206 is both a "composite" SNR, containing a 
PWN and shell structure together, and a "mixed morphology" SNR, with centrally concen- 
trated X-ray emission that cannot be simply associated with the PWN. Indeed, N206 is not 
the only SNR to show this combination; for instance, the Galactic remnant W28 combines 
centrally concentrated X-ray emission with a probable off-center PWN (Rho & Borkowski 
2002; Kassim 1992). Other examples include W44 (Cox et al. 1999; Rho et al. 1994; Harms, 
Hughes, & Helfand 1996) and IC443 (Kawasaki et al. 2002; Bocchino & Bykov 2001; Olbert 
et al. 2001; Petre et al. 1988). We shall call this category "mixed-composite" SNRs, reflecting 
the combined action of separate physical processes influencing their emission. 
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Table 1. Optical Echelle Observations 



Date Line Exposure Position 



2000 Dec 6 Ha + N II 1200 s E-W through "linear feature" 
2000 Dec 6 Ha + N II 1200 s N-S through SNR center 

2004 Jan 14 Ha + N II 1200 s E-W through bright filament 
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Table 2. X-ray Regions 



Label Name ACIS cts MOS cts Hardness 



Joint Chandra ACIS / XMM MOS regions 



1 


Whole SNR 


19,950 


10,920 


-0.49 


2 


Outer Limb 


7,450 


4,130 


-0.45 


3 


Central 


6,350 


3,320 


-0.58 


4 


Wedge 


2,375 


1,400 


-0.46 


5 


Linear Feature 


1350 


430 


-0.24 


Small 


spatial regions (Chandra ACIS 


only) 




6 


Knot 


456 




0.26 


7 


22" W of knot 


301 




-0.35 


8 


44" W of knot 


264 




-0.43 


9 


Point Src 


94 




0.50 


10 


Bow Shock 


754 




-0.45 


11 


Trail 


196 




0.076 



Note. - - Background-subtracted source counts. "Hardness" refers to the hardness ratio 
derived from Chandra ACIS counts. The "Hard" band is defined as 1.0-8.0 keV, the "Soft" 
band as 0.3—1.0 keV; the hardness ratio is H— S/H+S. 
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Table 3. Fits to X-ray emission for N206 SNR 



Spectral Nn 
Model cm~ 2 


kT 
keV 


abundance 
frac. solar 


T 

cm~ 3 s 


norm 

cm -5 


x 2 

reducd 


dof 




Outer Limb Region (excludes "linear feature") 


TV IT ^ 1 A /~1TO 

Merged ACIS 
nshock 2 7+ 1 1 v 1f) 21 
mekal 3.6t^ x 10 21 


n 47+0-04 
u - 4 ' -0.11 

22+ - 01 

u - zz -0.03 


u - zo -0.08 

n i4 +0 09 


1 4 +0 ' 9 x 10 11 

—0.4 


7±5 x 10~ 4 
6±2 x 10~ 4 


1.38 
1.93 


154 
155 


EPIC MUSl&z 
nshock 7 +5 x If) 20 
mekal 2.3181 x 1021 


o ^n+ 06 

24+ - 02 
u - z ^-o.oi 


o.2oig:S 
o.o8i8:g 


S 4+H x ID 11 


o_i_9 y in-3 

0.1±0.02 


1.17 
1.48 


83 
84 


ACIS + EPIC M0S1&2 
pshock 2.2±g;t x 10 21 

mekal 3.li^ x 10 21 


n 44+ - 04 

u - 4 *-0.06 

23+ 01 

u - zo -0.01 


n 9 9+0.04 
u - zz -0.04 

12+ - 05 

u - lz -0.02 


2.118:5 x io 11 


7±3 x 10~ 4 
1.5±0.3 x 10~ 2 
6±2 x IO" 3 
0.12 ± 0.03 


1.54 
1.96 


241 
242 


Central Region (excludes 


"linear feature") 










Merged ACIS 














pshock 3.5 ±0.4 x 10 21 
mekal 3.9 ±0.3 x 10 21 


0.34 ± 0.08 
0.22 ± 0.08 


0.25 ± 0.05 
0.4 ± 0.2 


4 ±2 x 10 11 


1.4 ±0.9 x 10~ 3 
3 ±2 x 10~ 3 


3.09 
3.79 


91 
92 


EPIC MOS1&2 














pshock 2.0 ±0.4 x 10 21 
mekal 2.7 ±0.3 x 10 21 


0.38 ± 0.04 
0.22 ± 0.01 


0.14 ± 0.03 
0.11 ± 0.02 


2.7 ±0.9 x 10 11 


2 ±1 x 10~ 3 
1.0 ±0.4 x 10~ 2 


2.10 
2.98 


64 
65 


ACIS + EPIC MOS1&2 














pshock 2.9±0.3 x 10 21 


0.37± 0.02 


0.25± 0.03 


3.0±0.2 x 10 11 


9±4 x 10~ 4 


3.21 


159 



2.1±0.9 x 10~ 3 

mekal 3.6±0.2 x 10 21 0.22± 0.06 0.4± 0.1 ••• 3±1 x 10~ 4 3.97 160 

6±2 x 10~ 3 



Note. - 



All spectra cover the range between 0.3-10.0 keV. 
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Table 4. Fits to X-ray Emission for "Whole SNR" Region of the N206 SNR 



Parameter 




Model Fits 








NEI 


Combined NEI + powerlaw 


component 


vmekal 


vpshock 


vpshock 


powerlaw 


N H (cm~ 2 ) 


2.81°,;} x 10 21 


2.21°;} x 10 21 


2.21°;} x 10 21 


kT (kev) 


n o/^/f +0.002 


0-453l°;°° 5 


0.431°;°} 




o/o 


0.36l°;° 2 


n OO+0.02 

u - oz -0.01 


0.301°;°} 




Ne/Ne 


24+ - 03 


n OO+0.02 
u -^-0.01 


0.241°;°} 




Mg/Mg 


u -°°-0.05 


0.391°;°^ 


41 + °-° 4 




Si/Si© 


0.918:1 


0.281°;°^ 


3 +ai 

u -°-0 .2 




Fe/Fe 


0.181°;°} 


u - zu -0.02 


21 + °-° 2 

u - z,± -0.02 




r (cm -3 s) 




3.5l°; 2 x 10 11 


3.3l°; 2 x 10 11 




r 








9 9 +0.3 
z - z -0.4 


ACIS norm (cm -5 ) 


5.251°;°? x 10- 3 


1.75l°;° 3 x 10- 3 


1.76l°;° 3 x 10- 3 


2.3±g;| x 10~ 5 


MOS norm (cm" 5 ) 


1.051°;°* x 10- 2 


3.471°;°^ x KT 3 


3.591°;°^ x 10- 3 


71? x 10~ 6 


Xrcd 


2.73 


2.23 


2.15 




dof 


294 


294 


292 





Note. — Results in Tables 4-8 are from simultaneous fits to ACIS and EPIC MOS data. 
All spectra cover the range between 0.3-8.0 keV. 
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Table 5. Fits to X-ray Emission for "Outer Limb" Region of the N206 SNR 



Parameter 




Model Fits 






CIE 


NEI 


Combined NEI + powerlaw 


component 


vmekal 


vpshock 


vpshock 


powerlaw 


N u (cm -2 ) 


2.9±g;| x 10 21 


1.7Iai x 10 


i.7l8:l x io 21 


trT (\rc\T\ 

Kl \K& V ) 


n OQ7+0.03 
U.Z.O 1 _ .02 


4fi+ - 02 


4fi+ 01 




o/o 


0.2518;" 


22 +0 - 01 


n 22 +0 - 01 




Ne/Ne 


o.2oi8:£ 


n 2 c i +0 - 03 

u - zd -0.03 


95 +0 - 02 




Mg/Mg 




n 4Q+0.06 
u -^°-0.06 


4S+ - 07 




Si/Si 




2 +0 - 2 


2 +0 - 2 




Fe/Fe 


o.2ii8:i 


n 24+ - 02 


24+ - 02 




r (cm -3 s) 




2.218:2 x io 11 


2.218:2 x io 11 




r 








2l 3 


ACIS norm (cm -5 ) 


2.69l8:8e x 10- 3 


5.H8:} x io- 4 


5.0I8: 2 x io- 4 


2l 4 x 10- 6 


MOS norm (cm~ 5 ) 


5.818:? x io- 2 


i.nl8:S x io- 2 




n +2xl0- 6 


Xied 


1.71 


1.47 


1.48 




dof 


203 


203 


201 
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Table 6. Fits to X-ray Emission for "Central" Region of the N206 SNR 



Parameter 


CIE 


Model Fits 
NEI Combined NEI + powerlaw 


component 
N u (cm -2 ) 


vmekal 
2.8l°;{£ x 10 21 


vpshock 
3.018;} x 10 21 


vpshock powerlaw 
2.7±g;| x 10 21 


kT (keV) 


o-278±8:8Sg 


n QQc+0.002 
u.ood_ 0006 


n 07+0.02 
u -° ' -0.01 




o/o 


0.5618$ 


4fi+ - 02 


0.46lg;« 




Ne/Ne 


o.iel™ 7 


24 +0 - 03 


0.2718;" 




Mg/Mg 


o.9±8:i 


S+ - 2 


0-618:1 




Si/Si 


n 7 +0 - 2 

u -'-0.3 


3 +0 - 2 

U.O_ 2 






Fe/Fe 


n 22+ - 01 


21 +0 - 02 


n 00+0.02 

u - zo -0.02 




t (cm -3 s) 




9+} x 10 11 


5.51} x 10 11 




r 








311 


ACIS norm (cm -5 ) 


1.22t°;°3 x 10- 3 


1.13lg;« x 10- 3 


7.7l°; 2 x 10- 4 


2±g : l x 10- 6 


MOS norm (cm -5 ) 


2.72l°;°« x 10- 3 




i.7i8;l x 10- 3 


n +3xio- 6 

u -0 


Xred 

dof 


2.50 
150 


2.22 

150 


2.18 
148 
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Table 7. Fits to X-ray Emission for "Wedge" Region of the N206 SNR 



Parameter 




Model Fits 






CIE 


NEI 


Combined NEI + powerlaw 


component 


vmekal 


vpshock 


vpshock 


powerlaw 


N H (cm" 2 ) 


2.318J x 10 21 


2.2+°-} x 10 21 

— U. 1 


1.818:} x 10 21 


kT (keV) 


o 32 +0 - 01 

u - oz '-0.01 


0.9018:81 


fifi+ 03 

u - uu -0.02 




o/o 


n 31 +0 - 04 

u - oi -0.04 


0.30^ 


u - oy -0.03 




Ne/Ne 


n nq+o.02 


n 1 O+0.04 
u - iz -0.04 


o.i4i8:8^ 




Mg/Mg 


o.36±g:Sg 


2 +ai 


o.4i8:l 




Si/Si© 




0.218:? 


o.H8:? 




Fe/Fe 


o-06±8:S 


o.i9l8:81 


o.i9l8:8I 




t (cm~ 3 s) 




71} x 10 10 


2.ll8i x 10 11 




r 








1 fi +0 - 9 

l -°-0.9 


ACIS norm (cm -5 ) 


5.9±g;g x 10~ 4 


1.0518:81 x io- 4 


i.i4i8:8i x io- 4 


2l 2 x 10~ 6 


MOS norm (cm -5 ) 


4.3t° 2 2 x ID" 2 


7.5l8i x io- 3 


7.918:1 x io- 4 


2.3181 x 10~ 4 


Xred 


1.49 


1.07 


0.98 




dof 


102 


102 


100 
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Table 8. Fits to X-ray Emission for "Linear Feature" Region of the N206 SNR 



Parameter 




Model Fits 






CIE 


NEI 


Combined NEI + powerlaw 


powerlaw 


component 


vmekal 


vpshock 


vpshock powerlaw 


powerlaw 


N u (cm- 2 ) 


i.oig;| x io 21 


3.5l°; 4 x 10 21 


2.2183 x 10 21 


2.618;^ x 10 2 


kT (keV) 


2-0l°j 


o O+0.9 

Z.O_Q g 


o.4l8:? 




o/o 


2l 3 


u - iJ -0.08 


0.218;^ 




Ne/Ne 


2tl 


o.is±8:8S 


o.2i8:^ 




Mg/Mg 


H? 


0.318;? 


0.7181 




Si/Si 


o i +0 - 5 

u - i -o.i 


o.2l8:i 


0.4183 




Fe/Fe 


1 +0 
u - 1 -o.i 


11+ - 05 

u - 11 -0.05 


u - 1 -o.i 




r (cm -3 s) 




1±\ x 10 10 


5xl0 13a 




r 






9 9+O.3 
z - z -0.3 


9 O+0.2 


ACIS norm 


4.51°;^ x 10" 4 


4.2l8;s x 10~ 5 


31^ x 10~ 5 i.5l8;3 x 10~ 5 


2.4l8:e x 10- 


MOS norm 


i-il8:2 x io- 4 


1.018;? x io- 4 


2.018:3 x io~ 4 2.7l8;7 x 10- 5 


6l£ x 10~ 5 


Xred 


1.13 


1.04 


0.93 


1.08 


dof 


92 


92 


90 


92 



a Error estimation gave a range of 6xl0 12 to 5xl0 13 cm 3 s; the latter is the hard limit for 
this parameter. At this upper limit the plasma is considered to be in ionization equilibrium. 
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Table 9. Small-region fits (Chandra ACIS only) 



Region 


kT 


T 


r 


K-pshock 


K 

- Lv pow 


X 2 /dof 


Knot 

22" W of knot 
44" W of knot 
Point Source 


0.918:1 
2.018:? 
o.36i°:8^ 


4±| x 10 10 
31} x 10 10 
811 x !0 n 


1-718:1 

21} 
31} 

9 i +0.4 
z - 1 -0.5 


2.618:9 x io- 6 

71} x 10- 6 

2.418:4 x 10- 5 


71} x 1(T 6 

1.2I8:? x io- 6 

1.818J x 10-6 
2.4l8:s x 1(T 6 


0.89/18 
0.67/15 
0.78/15 
1.4/6 



Note. - - All spectra cover the range between 0.3-8.0 keV. For these fits, combined NEI 
plasma (pshock) and power-law (pow) models were used. N# was fixed at a value of 2.2 xlO 21 
cm~ 2 and abundances were set to 30% solar, consistent with similar fits for this area above. 
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Table 10. Derived Physical Parameters for SNR N206 



Region 



Cool shell 


density 


10±4 cm -3 


Cool shell 


mass 


160 - 2000 M 


Cool shell 


expansion velocity 


202 ±5 km s~ 3 


Cool shell 


kinetic energy 


1 - 8 xlO 50 erg 


Cool shell 


thermal pressure 


3±1 x 10- 11 dyne cm" 2 


Hot gas 


temperature 


0.4±0.1 keV 


Hot gas 


density 


0.24 /- 1 / 2 cm- 3 


Hot gas 


mass 


270 ±10 flU M & 


Hot gas 


thermal energy 


6±1 x 10 50 fUt erg 


Hot gas 


thermal pressure 


3.3 ±0.4 x 10- 10 fUt dyne cm" 2 
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Table 11. Linear Feature Thermal vs Nonthermal X-ray Emission (Chandra ACIS only) 



Region Power-Law Flux Thermal Flux Flux Ratio 

erg cm' 2 s _1 erg cm -2 s _1 pwrlaw/thrml 



By fitted spectral model components 



Knot 2.2±0.6 x 1(T 14 2±1 x 10~ 15 0.9/0.1 

22" W of knot 7±1 x 10~ 15 7.5±0.5 x 10~ 15 0.5/0.5 



44" W of knot 


2.7±0.4 x 10' 15 7.9±0.7 x 10~ 15 


0.3/0.7 


By energy range 


(Thermal: 0.3-1.0 keV; Powerlaw: 2.0-i 


1.0 keV) 


Knot 


2.3±0.2 x 10~ 14 4.3±0.3 x 10~ 15 


0.8/0.2 


22" W of knot 


8.3±0.1 x lO" 15 7.3±0.2 x 10~ 15 


0.5/0.5 


44" W of knot 


8±1 x 10~ 16 7.3±0.1 x 10~ 15 


0.1/0.9 
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H 1 1 1 1 1 1 1 1 1 1 h -H 1 1 1 1 1 1 1 1 r- 

5 h 32™20" ib" 32^0" 50 s 40^ 5* t 32 m 20 B 10" SE^O'SO" 40^ 



a (J2000) 

Fig. 1. — Smoothed images, over the same spatial region: (a) XMM-Newton EPIC-pn 0.2- 
10.0 keV (b) XMM-Newton EPIC-MOS 0.2-10.0 keV (c) Chandra ACIS 0.2-0.9 keV (d) 
Chandra ACIS 0.9-3.0 keV (e) Chandra ACIS 3.0-8.0 keV (f) ATCA radio 6 cm. XMM- 
Newton images are adaptavely smoothed with the same scaling map; Chandra images are 
smoothed with a Gausian (a — 3). 
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Fig. 2. — Adaptively smoothed Chandra ACIS 0.3-10.0 keV image (grayscale) with ATCA 
6 cm radio contours overlaid. 
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a (JSOOO) 

Fig. 3. — Combined HST WFPC2 mosaicked images of N206 in Ha. 
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« (J3000) 



Fig. 4.— HST WFPC2 mosaicked images of the southwest side of N206 in (a) Ha, (b) [S II], 
(c) [0 III] and (d) [S II] /Ha ratio map. 
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Fig. 5.— Chandra ACIS (crosses) and XMM-Newton EPIC-MOS (dotted crosses) spectra, 
fits and residuals for four extraction regions: (a) Region 2, outer limb; (b) Region 3, central; 
(c) Region 4, wedge; (d) Region 5, linear feature. 
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Fig. 6. — Smoothed Chandra ACIS image with regions for spectral analysis, as listed in 
Table 2, overlaid and labeled, (a) Region 1; (b) Region 2; (c) Regions 3-4; (d) Region 5; (e) 
Regions 6-9; (f) Regions 10-11. 
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Fig. 7. — Top: echelle spectrum along "linear feature" . Bottom: ATCA 6 cm radio image 
of "linear feature" showing slit location. 



